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Dileptons, as an electromagnetic probe, are crucial to study the properties of a Quark-Gluon 
Plasma (QGP) created in heavy ion collisions. We calculated the invariant mass spectra and the 
anisotropic emission of thermal dielectrons from Au+Au collisions at the Relativistic Heavy Ion 
Collider (RHIC) energy ^snn = 200 GeV based on EPOS3. This approach provides a realistic 
(3+l)-dimensional event-by-event viscous hydrodynamic description of the expanding hot and dense 
matter with a very particular initial condition, and a large set of hadron data and direct photons 
(besides V 2 and vs !) can be successfully reproduced. Thermal dilepton emission from both the 
QGP phase and the hadronic gas are considered, with the emission rates based on Lattice QCD and 
a vector meson model, respectively. We find that the computed invariant mass spectra (thermal 
contribution + STAR cocktail) can reproduce the measured ones from STAR at different centralities. 
Different compared to other model predictions, the obtained elliptic flow of thermal dileptons is larger 
than the STAR measurement referring to all dileptons. We observe a clear centrality dependence of 
thermal dilepton not only for elliptic flow V 2 but also for higher orders. At a given centrality, of 
thermal dileptons decreases monotonically with n for 2 < n < 5. 

PACS numbers: 


I. INTRODUCTION 

Dielectron spectra have been measured by the 
PHENIX Collaboration [Ij and the STAR Collaboration 
@1 in Au+Au collisions at ^Jsnn = 200 GeV at Rel¬ 
ativistic Heavy Ion Collider (RHIC). In fact, dileptons 
have been proposed as a kind of clear probe in the study 
of the hot and dense matter named Quark-Gluon Plasma 
(QGP), created in heavy ion collisions [^, due to the fact 
that dileptons do not interact strongly and therefore sur¬ 
vive the whole evolution of a heavy ion collisions. This 
concerns the initial hard scattering Drell-Yan pairs, the 
thermal emission from the QGP and the Hadron Gas 
phase, and the decay of light and heavy flavor hadrons. 
The dielectron spectra from both STAR and PHENIX 
show the low mass enhancement beyond the cocktail. 
This phenomenon was also observed in the CERES di¬ 
electron i,i and NA60 dimuon data Medium mod¬ 
ifications of the p meson spectral functions have been con¬ 
sidered to describe this enhancement, based on hadronic 
many body effective theories 01| and using a relatively 
simple expanding fireball model. Another approach used 
the forward scattering amplitude method [T^, [HI and 
viscous hydrodynamics with smooth initial conditions as 
medium. These two methods could get a similar result 
[T^ . The combination of these two methods [T^ could 
also describe the STAR dielectrons data reasonably well, 
using a 2+1 dimensional ideal hydrodynamic model. In 
addition, a microscopic transport model- Parton-Hadron 
String Dynamics (PHSD) |14l - llq has successfully de¬ 
scribed the observed dimuon enhancement by NA60 and 
dielectron enhancement by STAR. 

Elliptic flow is the second Fourier coefficient of the az¬ 


imuthal particle distribution. The thermal elliptic flows 
of dileptons has been studied in hydrodynamic simula¬ 
tions with smooth initial conditions and also with 

fluctuating initial conditions (T^ . The fluctuations seem 
to play an important role in elliptic flows of dileptons in 
heavy ion collisions. However, the predicted elliptic flows 
of thermal dileptons are small, so do their predicted el¬ 
liptic flows of direct photons. The observed large elliptic 
flow of direct photons seems a big puzzle in our field. Us¬ 
ing EPOS3102, a (3+I)D viscous hydrodynamic model 
well constrained by hadron data, it seems that one can 
reproduce naturally the measured elliptic and trianger 
flow of direct photons (^. In this paper we will compute 
thermal dilepton production in an approach which has 
been successfully tested concerning hadron and photon 
production. This will serve as a benchmark to a system¬ 
atic compilation of hadron, dilepton and photon results. 


In additional to the thermal sources, dileptons may 
be produced from non-thermal sources such as primor¬ 
dial Drell-Yan annihilation and electromagnetic final- 
state decays of long-lived hadrons. To account for these, 
STAR cocktail results from Refs Q have been employed 
in this paper. 


The paper is organized as following. In section 2, we 
briefly review the calculation approach of thermal dilep¬ 
ton emission rates from quark-gluon plasma (QGP) and 
hadronic matter. In section 3, the results will be pre¬ 
sented. Finally, a conclusion is made in section 4. 
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II. CALCULATION APPROACH 

The calculation of dilepton production is obtained as 
an integration over an emission rate ^ = d^xdAq ' based 
on the space-time evolution of a QGP and the hadronic 
stage according to EPOSSjllj. We use a QGP emis¬ 
sion rate from extrapolating lattice QCD results to h- 
nite momenta. The calculation of the hadronic emis¬ 
sion rate is done by employing the vector meson domi¬ 
nance model. We consider explicitly the p meson, for the 
rest we refer to the so-called cocktail according to the 
STAR experiment (dileptons from decays of measured 
resonances) This procedure to compute dileptions, 
discussed in more detail in the following, will be referred 
to as EPOS3+DiLl. 


temperature T(x) of matter for a given space-time x, 
starting from some initial proper time tq. This is the 
basis of the dilepton (rate) calculations discussed later. 
In the hadronic phase, we have from EPOS a complete 
description of hadron trajectories. We may use this to 
compute energy densities, flow velocities and net baryon 
density. Assuming that the system may be approximated 
by a resonance gas in equilibrium, we use the correspond¬ 
ing equation of state, c.f. Apendix G of [l^, to get the 
temperature T and baryon chemical potential /is at any 
point in space time. This will be used for the dilepton 
rate calculations. 


B. Invariant mass spectrum 


A. The EPOS3 dynamics 

As explained in (2l|, EPOS3 is an event generator 
based on a 3+lD viscous hydrodynamical evolution start¬ 
ing from flux tube initial conditions, which are generated 
in the Gribov-Regge multiple scattering framework [^ . 
An individual scattering is referred to as Pomeron, iden- 
tihed with a parton ladder, eventually showing up as flux 
tubes (or strings). Each parton ladder is composed of a 
pQGD hard process, plus initial and final state linear 
parton emission. Nonlinear effects are considered by us¬ 
ing saturation scales Qs, depending on the energy and 
the number of participants connected to the Pomeron in 
question. 

The final state partonic system (corresponding to a 
Pomeron) amounts to (usually two) color flux tubes, be¬ 
ing mainly longitudinal, with transversely moving pieces 
carrying the pt of the partons from hard scatterings 
1^ . One has two flux tubes, based on the cylindri¬ 
cal topology of the Pomerons, but each quark-antiquark 
pair in the parton ladder will cut a string into two, in 
this sense one may have more than two flux tubes. In 
any case, these flux tubes constitute eventually both bulk 
matter, also referred to as “core” [ 2 ^ (which thermalizes, 
flows, and finally hadronizes) and jets (also referred to as 
“corona”), according to some criteria based on the energy 
of the string segments and the local string density. 

Concerning the core, we use a 3-l-lD viscous hydrody¬ 
namic approach, employing a realistic equation of state, 
compatible with IQCD results. We employ for all calcu¬ 
lations in this paper a value oi rj/S = 0.08. Whenever a 
temperature of Th = 168 MeV is reached, we apply the 
usual Cooper-Frye procedure to convert the fluid into 
particles. From this point on, we apply a hadronic cas¬ 
cade ( 2 ^, based on hadronic cross sections. 

So we have two phases, a plasma phase with a hydrody¬ 
namical expansion, and a hadronic phase, with individual 
hadron-hadron collisions. From the fluid dynamical ex¬ 
pansion of the plasma, we get the complete space-time 
information, i.e. the collective velocity v(x) (which de¬ 
fine the local rest frame (LRF) at each point x) and the 


The invariant mass spectrum of thermal dileptons is 
calculated as an integration of the dilepton emission rates 
over the realistic space-time evolution of the heavy-ion 
collision system and the momentum space of the dilep¬ 
tons 0 , namely. 



dNu 

d*xd*q 


{T{x),q*). 


( 1 ) 


Here, M and q are the invariant mass and four momen¬ 
tum of the dilepton pair, satisfying = q"^. The four 
momentum q* of the dilepton pair as viewed from local 
rest frame is given as 


Qo = 790 - 79 • d 


( 2 ) 


and 


q* =q+{l-l){q-v)-v/v‘^-jq°v, (3) 

where v is the flow velocity, = v ■ v and 7 = 
l/-\/l — jcd-. The flow velocity v{x) and the temper¬ 
ature Tix) at each space-time point x are provided by 
the above-mentioned EPOS3 dynamics. The momentum 
integral ranges from 0 to 5GeV/c (results do not depend 
on the upper bound as long as it is much bigger than 
the temperature). The quantity d^x = Tdrdxdydrj is the 
space-time element. 

The space-time integral starts from the initial time tq 
of the hydrodynamical evolution up to the end of hadron 
cascade. In the hadronic phase, uRQMD has been em¬ 
ployed after freeze-out in EPOS, and macroscopic quan¬ 
tities such as energy-momentum and flow velocity are 
reconstructed in order to calculate thermal emissions of 
dileptons and photons, as mentioned earlier. The emis¬ 
sion rate is the number of dileptons emitted per space- 
time unit and per energy-momentum unit. Thermal 
dileptons emitted from both the QGP phase and the 
hadron gas phase are considered, as explained in the fol¬ 
lowing. 
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C. Dilepton Emission Rate in the Quark-Gluon 
Plasma 


For dileptons emitted from the QGP phase, asymp¬ 
totic freedom implies that the production rate in the in¬ 
termediate invariant mass region (1. < M < 3.0GeP/c^) 
at high temperatures and densities can be described by 
perturbation theory as in Ref.Q. Recent progress using 
thermal lattice QCD (IQCD) to calculate dilepton rates 
nonperturbatively at vanishing three-momentum q = 0 
has been reported in Ref. [^. For practical applica¬ 
tions, an extrapolation to finite q values is needed. A 
construction proposed in Ref. Q and employed in our 
calculation reads 


Is 


= ^f^iqo, r)GEM{l + — ln[^^] 

47r^ q 1 + X- 


T'2 

+ 2 ™, —iFR(M 2 )ln(l 


2.912go 2gg+M% 
4,7: a sT Zq^ ’ 


(4) 


where a is the electromagnetic coupling constant, 
f^{qQ,T) the thermal Bose distribution. Gem = 
J2q=u,d,s 2 ;± = exp[-{qQ±q)/2T]. The quantity as is 
the temperature-dependent strong coupling constant, K 
a constant factor (equal 2), and finally we have a form 
factor with A = 2T. 


D. Rates from hadronic medium 


According to the vector meson dominance model |^ . 
the hadronic electromagnetic current operator is equal to 
the linear combination of the known neutral vector meson 
field operators, most notably V=^ w, 4>. This describes 
dilepton production successfully [3|, where the dilepton 
emission rate via the vector meson V is na 


SRy 

Sq 


a my 


■ItciD 


1 


V ^ 


-iii 1 ' 

e T — 1 


(5) 


The vacuum part TV’S is obtained from Gounaris- 
Sakurai formula as Refs [2^, [s^ , which could describle 
the pion electromagnetic form factor well, as measured 
in e+e“ annihilation Refs [sij: 


Ren™= = 


gj^ 

487r^ 




1 + t/l- 


4m2 

IvP" 


1- a/1- 


4m^ 

~JP~ 


( 8 ) 


IQ 2 ( ^ ^ w'^0 +P0X1 

M ^ Wo 


ImR^^ = - 




487r 


.(1 _ AZZL) 
^ M 2 ’ 


(9) 


Here, 2wo = Wp = 2 Srn^ + Pq■ vacuum width 


IS 


T^vac _ _i£_m ( P2.\3 
^ p ~ 4871- "‘‘P\ujn > ■ 


The interactive part Iipa{E,p) is obtained from p scat¬ 
tering from hadron of type a in the hadronic gas 


n. 


iiE,p) = -47r J 


dS 

(27r)- 


( 10 ) 


w 


where is the forward scattering amplitude in the 

c.m. system, E and p are the energy and momentum of 
the p meson, w2 = nS + S. The most copious hadrons in 
the hadronic gas such as tt mesons (a = tt) and nucleons 
(a = JV) are considered in our calculation. The quantity 
Ua is the Bose-Einstein occupation number of tt mesons 
and the Fermi-Dirac occupation number for nucleons, 
with temperature T and baryon chemical potential ps of 
the thermal bath provided by EPOS3 mentioned above. 
According to Eq.(I2]131), the emitted dileptons will get a 
Lorentz boost with the collective flow velocity offered by 
EPOS3, so the interactive term of self-energy and the re¬ 
sulted emission rate are calculated in the rest frame of 
the thermal bath. 

The pa forward scattering amplitude in the center-of- 
mass frame could be written as 


with the coupling constant gv determined by the mea¬ 
sured decay rate of vector-meson to dilepton production 
in vacuum, and my being the mass of the vector-meson. 
The imaginary part of the retarded vector meson propa¬ 
gator is given as 

^ (M2 -ml- Renv)2 + (Imnv)2 ’ 

where Hy is the self-energy of the vector meson V. 

We consider V=p and neglect the complexity of w and 
which seems consistent with the STAR dilepton data 
taking j2|. The p meson self-energy 

np = n;- + ^np„, (7) 

a 

for the contribution from the vacuum and from p meson 
scattering with hadrons of type a in the hadronic gas, 
respectively. 


/pa™' = (1 - a;)/Res + a;/Reg + /pom, (H) 

where the Pomeron term is dual to the background upon 
which the resonances are superimposed, cc is a function 
that matches the low energy Breit-Wigner resonances 
and high energy Reggons (dual to s-channel resonances) 
smoothly: 

a; = 0.5(1-f tanh (;^^—^——)), (12) 

where Aa is IGeV and 4GeV for pi: and pN scattering, 
respectivelvfl^. 

Regge and Pomeron term have the same form [ 13 : 


/r. 


,eg/P om 


4773 sin TTa ’ 


(13) 


where the intercept a and residue are 0.642, 

28.59, 12.74 for Regges and 1.093, 11.88, 7.508 for 
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Pomerons, respectively (the units yield a cross section 
in mb with energy in GeV) . 

The resonance term is |l0| : 


fues — 


‘^Qc.r 


E yt^R _ ^ R^pa _ 


(14) 


which involves a sum over a series of Breit-Wigner reso¬ 
nances of mass Mu and total width T/j. The resonances 
R used in our calculation are listed in table I for a = tt 
and table II for a = N, with R’s name, mass, decay 
width, branching ratio, spin, isospin and relative angular 
momentum Ref. (l^ . . 


Table I: Meson resonances R for pvr processes. 


Nanie(R) 

Mass(Mi{) 

r 

BR 

S 

IS 

L 

0(1020) 

1.020 

0.0045 

0.13 

1 

0 

1 

hi{mo) 

1.170 

0.36 

1 

1 

0 

0 

ai(1260) 

1.230 

0.40 

0.68 

1 

1 

0 

7r(1300) 

1.300 

0.40 

0.32 

0 

1 

1 

02(1320) 

1.318 

0.107 

0.70 

2 

1 

2 

a;(1420) 

1.419 

0.174 

1 

1 

0 

1 



Table II: Baryon resonances R for pN processes. 


Name(R) 

Mass(Mii) 

r 

BR 

S 

IS 

L 

N(2090) 

1.928 

0.414 

0.49 

0.5 

0.5 

0 

N(1700) 

1.737 

0.249 

0.13 

1.5 

0.5 

0 

N(2080) 

1.804 

0.447 

0.26 

1.5 

0.5 

0 

N(2190) 

2.127 

0.547 

0.29 

3.5 

0.5 

2 

N(2100) 

1.885 

0.113 

0.27 

0.5 

0.5 

1 

N(1720) 

1.717 

0.383 

0.87 

1.5 

0.5 

1 

N(1900) 

1.879 

0.498 

0.44 

1.5 

0.5 

1 

N(2000) 

1.903 

0.494 

0.60 

2.5 

0.5 

1 

A(1900) 

1.920 

0.263 

0.38 

0.5 

1.5 

0 

A(1700) 

1.762 

0.599 

0.08 

1.5 

1.5 

0 

A(1940) 

2.057 

0.460 

0.35 

1.5 

1.5 

0 

A(2000) 

1.752 

0.251 

0.22 

2.5 

1.5 

1 

A(1905) 

1.881 

0.327 

0.86 

2.5 

1.5 

1 

N(1520) 

1.520 

0.124 

0.008 

1.5 

0.5 

0 

A(1232) 

1.232 

0.118 

0.006 

1.5 

1.5 

1 


The c.m. amplitude and the scattering amplitude 
in the rest frame of a, /pa, are related by 

V~sfpa"'is)=mafpa{Ep), (15) 

with s = m^ + m^ + 2Epma. The latter is plotted in FigH] 
and FigUl for pir and pN scattering, respectively. The 
total contribution (Resonances + Pomeron background + 
Regge) in both are shown as black solid lines, and both 
are consistent with Ref. M- Each individual process 
pTT —>■ i? is also shown in FiglT] Their sum is shown as 


Figure 1: (Color Online) The individual and total amplitude 
for pvr —>■ i? scattering (imaginary part: upper panel, real 
part: lower panel). 

black dashed line. We can see that the contribution of 
resonances dominants the low energy scattering, and the 
dominant channels are R = ai(1260) and hi (1170) due to 
their large branch ratios. In Fig[21 the relevant individual 
channels pN —>■ R are shown. The black dashed line is 
the total baryon resonances contribution, among which 
R = A(1905), A^(2000) and A(1720) dominate both the 
imaginary and the real part of the the amplitude for pN 
scattering. 

Putting the forward scattering amplitutes together, 
we get the interactive term of self-energy of p-meson. 
The imgaginary and real parts are shown in FigI31 for 
a fixed p-meson three-momentum q = 0.3 GeV/c at 
T = 150 MeV, with pir scattering dashed lines and pN 
scattering dotted lines, respectively. The contribution 
from vacuum is plotted as solid lines. The total one is 
plotted as dashed-dotted lines. The vacuum contributes 
when M > 2m^, then dominates over a wide region. 
The in-medium effect from interactive term seems rather 
small. 

Yet, the interaction with the medium gives a broaden¬ 
ing effect to the spectral density, p{q) = —lmDp{q)/TT. In 
FiglU ImUp is plotted as a function of invariant mass M 
for a p-meson momentum oi q = 0.3 GeV/c, in a thermal 
hadronic gas of temperature T=100, 150 and 200MeV, re¬ 
spectively. The vacuum contribution by itself is plotted 
as a red solid line. With the increase of medium tempera- 












































5 



Figure 2: (Color Online) The relevant individual and total 
amplitude for pN —^ R scattering (imaginary part: upper 
panel, real part: lower panel) . 




M (GeV) 


Figure 3: The self-energy of p-meson . 



Figure 4: (Color Online) The imaginary part of the p meson 
propagators as a function of invariant mass for a momentum 
of 300MeC/c and a nucleon density of njv = 0. Results are 
shown for vacuum(red solid line) and three temperatures, T = 
lOOMeC/c (black solid line), T — IbOMeVjc (black dashed 
line), T = 2QQMeVjc (black dashed-dotted line). 


ture, the p-mesoii scattering with hadron a in the medium 
broadens the spectral density more and more, but the 
peak remains close to p-meson mass. (Note: T=200MeV 
is only to visualize the broadening effect, though it is too 
high for hadronic gas. ) 

The emission rates at different phases are now obtained 
via Eq.(|3][5]). In FigO the emission rate from the QGP 
phase is plotted as red thin lines (dotted: 200MeV, solid: 
150MeV), from the HG phase as green thick lines (dotted: 
150MeV, solid: lOOMeV). A higher temperature makes a 
stronger di-electron emission. Nevertheless the HG phase 
provides a pronounced peak around the p-meson mass 
775MeV which exceeds the QGP contribution. This peak 
is essentially due to the vacuum term in the self-energy 
(black dashed line). The emission rate from HG phase at 
150MeV obtained from effective interaction Lagrangians 
( 3 ^ is also shown as a blue dotted dashed line. It is inter¬ 
esting to see the two rates are quite close to each other, 
though different approaches and channels are considered. 

We should mention here that the above-mentioned 
emission rates work typically for ideal hydrodynamics. 
For viscous hydrodynamics such as FPOS3, a viscous 
correction is needed. The effect of shear viscosity to the 
spectra and elliptic flow of dileptons in QGP phase and 
HG phase has been investigated]^. Similar work has 
been done for direct photons [^. The elliptic flow of 
dileptons seems more sensitive to viscous effect than their 
spectra, however, it remains a modest effect. 


E. Elliptic flow V2 and Triangular flow U3 

The elliptic flow V 2 and the triangular one V 3 of ther¬ 
mal dielectrons are calculated in a similar way as in Refs 
(35| . The azimuthal angle dependence of the invariant 
mass spectrum of dielectrons for a given event can be 
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Figure 5: (Color Online) Emission rates of dielectrons from 
QGP phase (red thin lines) at temperature: ISOMeV (solid 
line), 200MeV (dashed line), and HG phase (green thick lines) 
at temperature: lOOMeV (solid line), ISOMeV (dashed line). 
The vacuum contribution to HG rate at ISOMeV is shown 
as black dashed line. HG rate from effective interaction La- 
grangians at ISOMeV is shown as blue dashed-dotted line. 


decomposed into harmonics of the azimuthal angle (j) as 

^ ^ ■;^[l-|-2ti2COs2(^-V'2) + 2u3COs3((/>-'03)-|-], 

d(f> ZTT 

( 16 ) 

where V 2 {vn) is the elliptic flow (higher order harmonics), 
and 'ipn is the nth-order event plane angle. Obviously, Vn 
and ij^n depend on the dielectron’s invariant mass M and 
vary event by event. From Eq. dm), one can easily get 


, . , /o"cos(n(^)^|ed</) 

Ve,nCOS{n'tpe,n) = 


fc 


2-n- ^1 
0 d(/> 


dcj) 


(17) 




Ue,nSin(n^/>e,„) = 


(18) 


with the subscript e added to emphasis variables for a 
single event. Let’s note their right sides as < cos ncj) >e 
and < sinncj) >e, respectively. Then, for each event, the 
harmonics and reaction plane of order n can be obtained 
as 


and 


Ve,n = \/< COs{n(j)) >e + < sm{n(j)) (19) 

< sm{n(j>) >e 


'06,n = — arctan 


n < cos(n0) >e ' 


( 20 ) 


The elliptic flow V 2 and higher order harmonics of the 
sample is obtained via event average 


N 


E oe,n 


This is equivalent to the definition in 


( 21 ) 


III. RESULTS 

The results of thermal dielectrons calculated based on 
EPOS with initial time tq = 0.35/m/c will be shown in 
the following. 

In Figini is shown the invariant mass spectrum (left 
panel) and elliptic flow (right panel) of thermal di¬ 
electrons at rapidity y = 0 from AuAu collisions at 
y/SNN = 200 GeV with centrality 10-40%. The thermal 
results (black solid lines) have been decomposed into two 
contributions, coming from the QGP phase (red dotted 
lines) and the hadronic phase (blue dashed lines). 

From the left plot of Figl^l it is clear that the QGP con¬ 
tribution dominates the thermal spectra at M>lGeV/c^ 
region, as in the above figure on emission rates. And the 
peak from the contribution of the hadronic phase indeed 
remains around the p mass after the space-time integral 
of the emission rate. Thus, in the total thermal spectrum, 
a pronounced peak still exists. 

On the right panel of FiglHl the elliptic flow of dilep¬ 
tons from the QGP phase (dotted line) and from the HG 
(dashed line) are shown. The elliptic flow of the latter 
phase is much larger due to a longer development of ra¬ 
dial flow. The result of thermal dileptons (solid line), 
the average of the two flows weighted by their yields in 
each phase, ranges between the upper bound (values in 
hadronic phase) and lower bound (values in QGP phase), 
and is close to the corresponding bound of the dominat¬ 
ing phase. Therefore, no evident peak for the elliptic 
flows of dileptons from either phase still produces a peak 
of thermal dileptons around the p-meson mass. 

A similar phase competition in the invariant mass spec¬ 
trum and V 2 occurs at other collision centralities, ie, 0- 
10% and 40-80%. 

After investigating the thermal contribution, we would 
like to compare with data. Thus a study of non-thermal 
contribution is needed. EPOS3 is almost ready for this 
goal. For this moment, the employment of STAR cocktail 
data can also make an independent check of the thermal 
emission. STAR cocktail contribution includes all the 
non-thermal contributions, such as contributions from 
Drell-Yan process and decay of long lived hadrons, such 
as light mesons, D-mesons and so on. Based on the mea¬ 
surement of those hadrons, Tsallis fitting and decay sim¬ 
ulation, STAR offers the invariant mass spectra of the 
cocktail contribution [2^. 

In Figl3 we show the EPOS3+DiLl thermal con¬ 
tribution (dotted line) and the STAR cocktail (dashed 
line) from 0-10% Au+Au central collisions at ^/snn = 
200 GeV. Their sum as the invariant mass spectrum of 
dielectrons (solid line) is also compared to STAR data [1| 
(empty dots). 

In Fig[5Kleft panel) the invariant mass spectra of di¬ 
electrons from Au+Au collisions at -^snn = 200 GeV for 
different centralities from EPOS3+DiLl are compared to 
STAR data iia, from top to bottom: 0-10% (central), 
10-40%, 40-80%, and 0-80% (MinBias). For better vis¬ 
ibility, the latter three results are multiplied by factors: 
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Figure 6: (Color Online) Thermal di-electrons (black solid lines) from AuAu collisions at at =200 GeV with centrality 
10-40% are decomposed into di-electrons from QGP phase (red dotted lines) and hadronic phase (blue dashed lines). Left 
panel: the invariant mass spectra. Right panel: the elliptic flow. All calculations are based on EPOS3+DiLl. 



Figure 7: (Color Online) Invariant mass spectra of dielec¬ 
trons from Au+Au collisions at ^/sNN = 200 GeV for cen¬ 
trality 0-10%. The STAR dataQ (empty dots) are compared 
to the total contribution (solid line), which is the sum of the 
EPOS3+DiLl thermal contribution (dotted line) and STAR 
cocktail (dashed line). 


0.05, 0.02, 0.0002, respectively. The right panel shows 
the ratios of measured data and our calculations to the 
cocktail at different centralities. This precise comparison 
shows a deviation from the experimental data in low mass 
region for the first centrality. The calculated results de¬ 
scribe the data well for centralities of 10-40% and 40-80% 
better than 0-10%. 

In Figini the red solid line is our calculated elliptic 
flow of thermal di-electrons from Au+Au minimum bias 
collisions (0-80% centrality) at ^/snn = 200 GeV. As a 
reference, we also show STAR data (stars) of the 
elliptic flow of di-electrons (from all sources, of course). 
The calculated elliptic flow of thermal dielectrons is com¬ 


parable, in fact even larger than data for all di-electrons. 
The cocktail contribution will inherit a certain elliptic 
flow because hadrons carry elliptic flow before decaying 
into di-electrons. Most models predict it to be lower than 
the STAR reference [^, [T^, [I^ , whereas our result based 
on EPOS3 is larger. 

To make a complete collection of our theoretical re¬ 
sults, we present the centrality dependence and the n 
dependence of the flow harmonics Vn of dileptons in the 
following. 

In Fig[T0] is shown the elliptic flow of thermal dielec¬ 
trons from AuAu collisions at y/SNN = 200 GeV for dif¬ 
ferent centralities: 0-10%(pink dashed line), 10-40%(red 
dotted line), 40-80%(dark blue dashed-dotted line). The 
centrality dependence of high order n is also investigated. 
A strong centrality dependence, from central to periph¬ 
eral collisions, occurs not only to the elliptic flow but also 
to higher order, n = 3,4, 5. 

In FigfTTl the higher order harmonics coefficients 
(n = 2, 3,4, 5) of thermal dielectrons from An-An colli¬ 
sions at s/snn = 200 GeV with centrality of 10-40% is 
presented. The curves have same shape, same trend as 
direct photons. The magnitude decreases monotonically 
with the order n. And it is the case for all the investi¬ 
gated centrality class, 0-10%, 10-40% and 40-80%. The 
low pt charged hadrons behaves similarly, as observed by 
ALICE dil. 


It’s often said that the development of anisotropic flow 
is controlled by the anisotropies in the pressure gradients 
which in turn depend on the shape and structure of the 
initial density profile. The latter can be characterized by 
a set of harmonic eccentricity coefficients e„ as defined 
in (^ . The initial space eccentricity e„ as a function of 
centrality for n = 2, 3,4, 5 in EPOS3 is shown in FigfT^ 
which looks quite similar to the initial space eccentricity 
in IP-Glasma, MC-Glauber and MC-KLN as shown in 

la. 














Figure 8: (Color Online) Left panel: Invariant mass spectra of dielectrons from EPOS3+DiLl compared to STAR data ilia, 
for centrality 0-10% (central), 10-40%, 40-80%, 0-80% (MinBias) (from top to bottom). Right panel: Ratios of STAR dilepton 
invariant mass spectra and our calculations to cocktail, for centrality 0-10% (central), 10-40%, 40-80%, 0-80% (MinBias) (from 
top to bottom). 



Figure 9: (Color Online) Thermal elliptic flow of dielectrons 
from ^^snn = 200 GeV Au+Au minimum bias collisions and 
comparison with STAR data 


What builds up such a large elliptic flow of thermal 
dileptons, as shown in Fig. inp The same question is 
asked in the direct photons paper based on EPOS3 (^ . 


As explained there, most models predict an elliptic flow 
too low, whereas EPOS3 can reproduce the measured 
elliptic flow. A systematic study [23| on the time evo¬ 
lution of space eccentricity and momentum eccentricity 
of the hot dense matter in EPOS3 shows that a correct 
interplay between the momentum eccentricity and radial 
flow is important to account for the elliptic flow of both 
hadrons and electromagnetic probes, dileptons and direct 
photons. Compared to the surface emission of hadrons, 
the overall emission inside the hot dense matter makes 
direct photons and dileptons more sensitive to the radial 
flow. Thus a large elliptic flow of direct photons and 
dileptons are obtained at the same time, with a larger 
mean radial flow and smaller momentum eccentricity. It 
should be said that we use the default EPOS3, with all 
pamaters fixed from hadronic observables. 


IV. DISCUSSION AND CONCLUSION 

We investigated the anisotropic emission of thermal 
dileptons from Au+Au central collision at ^/snn = 
200 GeV at the RHIC, accompanying the study of the 
direct photons in where a possible solution to the 
direct photon puzzle is provided. 

With a detailed explanation of dilepton emission rate 
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Figure 10: (Color Online) The elliptic flow (left) and trianger flow (right) of thermal dielectrons from AuAu collisions at 
^/snn = 200 GeV. 0-10%(pink dashed line), 10-40%(red dotted line), 40-80%(dark blue dashed-dotted line). 



Figure 11: (Color Online) Predicated harmonics coefficients 
Vn {n = 2, 3, 4, 5) of thermal dielectrons are shown by various 
curves for centrality 10-40%. 



Figure 12: (Color Online) Eccentricity e„ as function of cen¬ 
trality for n = 2, 3, 4, 5 with EPOS3. 


and a relatively brief explanation of the space-time evolu¬ 
tion within EPOS3, we first calculated the thermal con¬ 
tribution (from QGP and p mesons) to the invariant mass 
spectra at different centralities, comparing with the cock¬ 
tail contribution. The thermal contribution does play 
an important role in the low invariant mass region, ie, 
M < IGeV/c^, which offers a useful window to study 
the hot and dense matter created in relativistic heavy ion 
collisions, especially the early stage of the hot and dense 
matter. Our calculated spectra, together with STAR 
cocktail, could reasonable well reproduce STAR’S mea¬ 
sured invariant mass spectra of di-electron at different 
centralities. 

Then, we investigated elliptic flow and higher order 
harmonics of thermal dileptons. The elliptic flow of ther¬ 
mal dileptons is predicted to be quite large, larger than 
in other models, and comparable to the dilepton data for 
centrality 0-80%. 

Finally, we made predictions of the of thermal dilep¬ 
tons for 2 < n < 5 at three centrality classes, 0-10%, 10- 
40%, and 40-80%. The Vn of thermal dileptons is found 
to decrease with n monotonically for the three centrality 
class . The same centrality dependence for 2 < n < 5 
is observed with EPOS3, namely, more central collisions 
provide smaller Vn of thermal dileptons. Further experi¬ 
mental measurements are highly awaited to test our so¬ 
lution to the puzzle of direct photon elliptic flow and 
improve our understanding of the created hot and dense 
matter. 
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